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In  many  queueing  systems  (e.g.^  a computer  facility)^  the  server 


is  subject  to  breakdown.  If  the  breakdowns  are  unpredictable  in  nature 


and  the  server  is  not  replaced  or  repaired  until  a breakdown  occurs  then 


the  facility  will  be  unable  to  provide  uninterrupted  service  to  its 


customers.  In  such  a case  it  is  important  to  understand  how  the  break' 


downs  will  effect  the  system's  level  of  performance 


Avi-Itzhak  and  Naor  [l]  and  Gaver  [5]  obtained  the  steady  state 


quantities  L and  W for  an 


queue  whose  unreliable  server 


alternates  between  operational  periods  of  exponential  duration  and  failed 


periods  of  arbitrary  random  duration.  Although  general  with  respect  to 


the  distributions  of  the  service  and  repair  times^  the  results  have  the 


disadvantage  of  requiring  a stationary  Poisson  arrival  process. 

Yechiali  and  Naor  [7]  and  Fond  and  Ross  [2]  both  considered  varia 


tions  of  a single  server  exponential  queueing  model  in  which  the  arrival 


with  the  length  of  time  the  system  operates  with  pair  (A  ^ 


exponentially  distributed  with  rate  cO;  i = 1,  2.  In  [7]^  the  steady 


state  quantity  L was  derived  for  a system  with  infinite  queue  capacity 


while  in  [2]  the  steady-state  proportion  of  customers  lost  was  obtained 


for  a system  in  which  any  arrival  finding  the  server  busy  is  lost.  The 
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special  case  1^2  = ® results  in  models  of  queueing  systems  subject  to 
breakdown  and  having  nonstationary  Poisson  arrivals;  however^  all  prob- 
ability distributions  are  restricted  to  be  exponential. 

The  model  developed  here  is  an  important  one  that  cannot  be 
obtained  as  a special  case  of  the  models  in  the  literature.  In  particular^ 
consider  a single  server  queueing  system  with  the  following  characteristics: 

(a)  The  system  alternates  between  two  states:  the  operational 
state  and  the  failed  state  (also  referred  to  as  the  repair  state). 

(b)  When  operational^  the  system  functions  as  an  K/E^I  queue; 
that  is^  customers  arrive  according  to  a Poisson  process  with  rate  A, 
and  service  is  according  to  an  Erlang  distribution  with  mean  l/jj  and 
shape  parameter  k. 

( c)  If  service  to  a customer  is  interrupted  by  a breakdown, 
resumption  takes  place  as  soon  as  the  repair  period  ends  with  no  loss  of 
service  involved. 

(d)  Although  no  service  takes  place  during  the  repair  period, 

customers  continue  to  arrive  according  to  a Poisson  process,  but  now 
having  a rate  of  instead  of  A. 

(e)  The  duration  of  operating  periods  is  exponential  with  mean 
I/cCt  and  the  duration  of  repair  periods  is  Erlang  with  mean  l/cp  and 
shape  parameter  m. 


This  model  is  not  a special  case  of  [7]  because  it  permits  both  the  service  ! 
and  repair  times  to  have  Erlang  distributions  rather  than  restricting  f 
them  to  be  exponential.  Furthermore,  the  model  is  not  a special  case  of  j 


if 


2 


[l]  and  [3]  because  the  Poisson  arrival  process  is  nonstationary  vdien 
^ A.  A nonstationary  arrival  process  is  useful  in  many  practical 
situations  where  the  customers  are  aware  when  the  server  is  inoperable. 

In  such  cases^  one  expects  to  find  < A or  even  Aj^  = 0. 

The  constant  c in  assumption  ( e)  controls  how  rapidly  the  system 
occilates  between  the  operational  and  failed  states.  Holding  Ct  and  p 
constant  while  increasing  c has  the  effect  of  keeping  the  steady  state 
probabilities  of  being  in  the  operational  and  failed  states  constant 
while  at  the  same  time  increasing  the  frequency  with  vrfjich  the  system 
changes  states. 

The  remainder  of  this  paper  is  organized  as  follows:  Section  1 
analyzes  the  queueing  system  described  by  assumptions  (a)  - because 

Little's  formula  is  not  valid  for  this  model  unless  A^  = \ the  steady- 
state  quantities  L and  W must  be  derived  separately.  Some  special 
and  extreme  cases  of  the  general  model  are  considered  in  Section  2. 
Section  5 not  only  shows  that  both  L and  W are  decreasing  and  convex 
functions  of  c but  also  investigates  the  behavior  of  the  system  as 
c -*m.  The  paper  concludes  in  Section  1+  with  an  application  to  a pro- 
duction-storage system. 
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1.  STEADY  STATES  RESULTS 


Preliminary  Analysis 


Throughout  this  section,  as  well  as  Sections  2 and  1+,  the  constant  c 
will  be  assumed  to  equal  1,  This  eliminates  the  need  to  carry  c along 
in  all  the  derivations  when  it  is  only  relevant  to  the  analysis  in  Section  3. 

The  "method  of  phases"  (cf.  [4,  p.  l68])  provides  a convenient 
means  of  obtaining  the  steady  state  results.  It  is  well-known  that  an 
Erlang  random  variable  with  mean  \/\x  and  shape  parameter  k is  equivalent 
to  the  sum  of  k independent,  exponentially  distributed  random  variables 
each  having  the  same  mean  l/k^i.  Hereafter,  both  the  Erlang  service  times 
and  Erlang  repair  times  of  the  model  will  be  viewed  as  consisting  of  a 
series  of  identical  and  independent,  exponentially  distributed  phases. 

The  system  can  now  be  analyzed  as  a continuous  time  Markov  process 
with  states  {(i,j)|i  = 0,1,  ...,  m and  j = 0,  1,  2,  ...)  where  i = 0 
denotes  the  system  is  operational,  1 < i < m denotes  the  number  of  phases 
remaining  in  the  repair  process  until  the  system  becomes  operational,  and 
j denotes  the  number  of  service  phases  in  the  system  ( the  s'jm  of  the 
number  of  phases  remaining  for  the  customer  in  service  and  k times  the 
number  of  customers  in  the  queue).  The  transition  probabilities  are 
stationary  and  satisfy  the  Kolmogorov  differential  equations.  Furthermore, 
the  steading  state  probabilities  exist,  are  independent  of  the 

initial  state,  and  satisfy  the  following  balance  equations: 
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(A^+m?) 

PmO 

^ a 

PqO 

( 

:j  = 0) 

(la) 

(Aj^-i-m?) 

^iO 

= m? 

'’i-(-l,0 

( 1 < i < m- 1)  1 

:j  = 0) 

(lb) 

(A+O) 

Pqo 

= m? 

^10 

+ k^x  1 

[j  = 0) 

( Ic) 

(Aj^-t-in?) 

^mj 

= a 

Poj 

+ H Vi-k 

(j  > 0) 

(Id) 

(A^-nr?) 

Pij 

= m? 

^i-4-1,  j 

*• ’'i  Pij-k 

(j  > 0^ 

( le) 

( A-fOf  k^^) 

Poj 

= m? 

Plj 

♦ ‘’o,j-k  " '‘“'’O.jtl 

(j  > 0) 

(If) 

where  a negative  subscript  in  ( Id)  - ( If)  indicates  the  term  is  zero. 

Figure  1 contains  the  portion  of  the  Markov  chain's  state  transition 
diagram  corresponding  to  states  (i,j)  with  0 < i < m and  j > k. 

It  is  clear  from  the  figure  that  equations  ( 1)  can  be  interpreted  as 
requiring  the  mean  transition  rates  into  and  out  of  a state  to  be  equal 
at  steady  state. 

Pv  ■ Poj  Pf  ‘ Cl  Co  "ij'  "w  '■f 

are  the  steady  state  probabilities  of  the  system  being  operational  and 
failed,  respectively.  On  considering  the  underlying  two  state  (operational 
and  failed)  stochastic  process,  it  is  immediate  that 


P„  = P/(Q^P)  , 


p = 0/(04-?) 


1 ^ 


' I 


i 


t. 


Let  the  average  arrival  rate  and  average  service  rate  in  steady  state 
be  denoted  by  A = 3*’^*  4 = I^Py»  furthermore^  let  r = A/h. 

It  will  be  demonstrated  shortly  that^  as  is  often  the  case,  A < p is 
a condition  for  steady  state.  Note  that  each  of  the  quantities  p^,  p^, 

A,  p,  and  r would  be  independent  of  c even  if  the  temporary  assumption 
c = 1 were  dropped. 


The  Generating  Function 

Generating- function  techniques  must  be  used  to  further  analyze 
the  model  as  there  is  no  way  of  solving  (1)  in  a recursive  manner  to 
obtain  closed- form  expressions  for  the  Define  the  generating 

functions 


G J z)  = ^^0  Pi  j 1^1  < ’•j  f = 0,1,2, ...,m 

=Co‘^i(^)  Izl<l.  (2) 


Multiplying  each  equation  of  the  sets  ((la),  (Id)),  {(lb),  (le)},  and 
((Ic),  (If))  by  z^  and  summing  over  all  j yields,  respectively, 

G^(z)  = [a/(A^+mP-A^z'')]  Gq(z)  , (3) 

G^(z)  = [mp/(Aj^+mP-A^z*^)  ] G^^j^(  z)  , (1<  i ^m-l)  , (‘^) 

k^l 

Gq(z)  = [m^z  Gj^(z)  - Iqi  Pqq(  1-z)  ]/ [( A-fOflqi)  z-Az  -Iqi]  . (5) 
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Equations  (3)  and  ( can  be  used  recursively  to  express  G^{z)  in  terms 
of  Gq(z)  as 

Gj^(z)  = (Q/itiP)[in(3/(>vj^+inP-A^z‘^)]"*  Gq(z)  (6) 

and  (5)  can  be  rearranged  as 

G^(z)  = (it^z)'^  ([(A40!+kn)z  - - k|i]  Gq(z)  + kupQQ(l-z)).  (7) 

Equating  the  expressions  for  Gj^(  z)  in  (6)  and  (7)  and  solving  for 
Gq( z)  gives 

Gq(z)  = {ku  Poo(l-2)  [f(z)  ]"’]/D(z),  (8) 

where 

f(z)  = 1 + (Aj^/mp)  (1-z*^)  , (9) 

and  the  denominator  is 

D(  z)  = az  + [kpi  + Az*^'*'^  - ( A+Ctflqj,)  z ] [ f ( z)  ]”*  . (10) 

Using  (5)  and  (4)  recursively  results  in 

G^(z)  = {(o/n^)  ku  Poo(l-z)  [f(z)]^'^)/D(z)  ( 1 < i < m)  . (11) 

G(  z)  can  now  be  calculated  from  (2),  (8),  and  (11)  as 


I 
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I 

i 

I 


7 


i 


G(z)  = Pqq[n(z)/d(z)] 


(12) 


where 

N(z)  = Ml-z)  {[f(z)]”  ^ (a/me)  [f(z)]^)  . (15) 

Using  G(  1)  = 1 and  computing  lim^  ^ G(  z)  by  applying  L'Hopital's 
rule  to  (12)  yields,  after  algebraic  simplification. 

Substituting  ( li<-)  into  (12)  results  in  a final  expression  for  G(  z) , 
that  is 


G(z)  = (1-r)  [p/(a^P)]  [N(z)/D(z)]  . (15) 

Since  p^^  > 0,  (14)  also  verifies  the  previously  mentioned  condition 
for  steady  state,  A < p.. 

The  busy  fraction  p = 1 - P-o  equals  1 - G( O)  and  can  be 

evaluated  from  (15),  (9),  (10),  and  (15).  Provided  A^  > 0, 

p = 1 . (1-r)  [p/(afP)]  (1  + (a/A^)  [1  - (mp/(Aj^+mp))"’])  . 

/N  /V  ^ /N. 

In  general,  then,  p /t  (A/p)j  however,  if  Aj^  = 0,  p = (A/p)  does'  hold. 
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Recursions  for  the  (p..) 

Unfortunately^  no  simple  relationship  exists  relating  the 

to  p , p,- p However  the  (p..)  can  be  computed  efficiently 

00^  10^  ^ nO  ' ij 

from  (14)  and  the  set  of  recursive  equations  (for  j = 0,  1,  2, 

Pmj  = (^1+"^)'^  ^ ^iPmJ-k) 

Pij  = (A^+mP)'^  (n#  p^^j^  j PiJ-k^»  ( i=m- l,m-2, . . . , 1}  ( l6b) 

P0,j  + L = 4j+l-k  Pon  ^ “ ^=0  Pon  ’ ^ ^0  Pfn^  ^ 

where  a term  is  zero  if  it  has  a negative  subscript  and  the  lower  limit 
of  summation  in  the  first  term  of  ( l6c)  is  reset  to  0 if  it  is  negative. 
Equations  ( l6a)  and  ( l6b)  are  obviously  equivalent  to  (Id)  and  ( le) 
while  ( l6c)  follows  from  (If)  and  a simple  inductive  argument.  Of  course, 
the  steady- state  probability  of  having  n customers  in  the  system  is  given 
by  1-p  for  n=0  and,  for  n > 0,  by  ^>(  n- 1)  k+1  Pfj  ’ 


Computation  of  L 

The  computation  of  L and  L^,  the  steady- state  average  number 

of  customers  in  the  system  and  in  the  queue,  respectively,  require  some 

P P 

preliminary  results.  Let  L and  be  steady- state  notation  for  the 

average  number  of  customer  service  phases  in  the  system  and  in  the  queue, 

p 

respectively.  Furthermore,  let  L denote  the  average  number  of  service 


I 

I 

i 


.1. 


phases  remaining  for  the  customer  (if  any)  in  service.  The  following 
relationships  clearly  hold: 


= kL  , 

q q ’ 


P P P 
L = L + L 
q s 


L = L + p 

q 


(17) 

(18) 

(19) 


Substituting  ( I7)  into  (18)  and  solving  for  results  in 


= (1/k)  (L^  - L^)  . 


(20) 


Relationships  ( I9)  and  (20)  then  yield 


L = (1/k)  (L^  - L^)  + p . 


(21) 


P P 

The  problem^  then^  is  to  compute  L and  L^. 

p 

Now  L = G '(!),•  however^  evaluating  G'(l)  is  not  easy. 
Expressing  G( z)  = p^^  N(z)/d(z)  and  using  L'Hopital's  rule  twice  gives 


G’(l)  = Poo[N’’(l)*D'(l)  - N'(1).D"(1)]/2[D’(1)]  . (22) 


The  algebraic  manipulations  required  by  (22)  are  straightforward  but 
quite  long.  Because  they  would  occupy  several  pages,  the  computations 
are  omitted^  however,  they  result  in 
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L^  = 


l-r 


k+1  m+I 
2 2m 


kQ!A^(li-A+Aj^) 


(25) 


Obtaining  L requires  the  development  of  another  generating 
s 


function.  In  the  M/E^/l  queue  not  subject  to  breakdown,  given  that 
a customer  is  in  service  in  steady- state,  the  number  of  phases  remaining 
until  his  service  is  complete  is  equally  likely  to  be  1,  2,  ...,  or 
k (cf.  [k,  p.  169]).  However,  this  is  not  the  case  when  the  queue 
is  subject  to  breakdown.  Define  the  generating  functions 


Hi(y)  . Pjo  + 7^1  (5^,0  ’’i.jk.n'y"  - lyl  < 1,  i . o,i. 


H(y)  . »i(y)  , 


ly|  < 1 


(24) 


Clearly,  = H'(l).  As  demonstrated  in  the  Appendix,  the  lengthy 
derivation  of  H(y)  results  in 


H(y)  = (l-p)  + (p-r)y‘'  + rt[y(l.y'')]/[k(l-y)]]  . (25) 


Application  of  L'Hopital's  rule  twice  to  (25)  yields 


= H'(l)  = (p-r)k  + r[(k+l)/2]  . 


(26) 


Finally,  combining  (21),  (23),  and  (26)  results  in 


L = 


l-r 


. k+1  m+1  QA^(h4-A^-A) 

( l-r)  + p,  + 2m  ’ .2 

^(a^-^) 


(21) 


l| 
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Computation  of  W 


Because  Little's  formula  is  not  valid  for  this  model  unless  = X, 
L cannot  be  used  to  immediately  obtain  the  steady- state  expected 
value  of  the  time  a customer  spends  both  in  the  queue  and  in  service. 


Instead^  let  W = where  the  (i  = 1,  2,  5)  are  the 

expected  values  in  steady  state  for 

1.  The  time  that  elapses  from  the  customer's  arrival  until  the 
queue  first  becomes  operational;  if  the  queue  is  operational  \dien  the 
customer  arrives^  this  time  equals  zero. 

2.  The  time  that  elapses  from  when  the  queue  first  becomes  opera- 
tional until  everyone  present  when  the  customer  arrived  has  been  served;  if 
the  arriving  customer  finds  the  system  empty,  this  time  equals  zero. 

5.  The  time  that  elapses  from  when  the  customer's  service  begins 
until  it  ends. 

Each  of  the  T^  ( i = 1,  2,  5)  will  now  be  calculated. 

Applying  L'Hopital's  rule  to  (11)  yields  G^(  1)  = Q/tm(Ct4-P)] 
for  i = 1,  2,  ...,  m so  that,  given  the  system  is  failed,  the  number 
of  exponential  repair  phases  remaining  until  it  becomes  operational  is 
equally  distributed  over  1,  2,  ...,  m.  Hence,  given  the  system  is 
failed,  the  expected  number  of  repair  phases  remaining  is  [(nH-l)/23. 

Since  each  repair  phase  has  an  expected  length  of  (1/tnP), 


= [(nH-l)/2np]  [cV(a^-P)  ] 


T 


The  times  Tg  and  can  both  be  divided  into  two  parts:  time 

w 

during  which  the  system  is  operational  ( i = 2^  5)  and  time  during  which 

the  system  is  failed  (T^j  i = 2^5).  Clearly,  Tg  = and  = l/n. 

w 

During  the  expected  operational  times  (i  = 2,  J),  the  expected  number 
w 

of  failures  is  Since  the  expected  repair  time  for  each  failure  is 


VP,  = (O^P)T^  for  i = 2,  5.  Hence, 


Tg  = [(afP)/p]  [L  /k4] 


(29) 


= [(afp)/pu] 


(50) 


Combining  (23),  (28),  (29)  and  (30)  results  in 


W = (A) 


• 1,  r X i , , k+1  TtH-l  a 

(ir^)  * 2ir  • ' ♦ 2W  • 51 


Ai(n+Ai-A)  - Am-v  1 


>1! 


(51) 


Note  from  (27)  and  (3I)  that  L = AW  does  not  hold  unless  A^^  = A. 
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2.  SPECIAL  AND  EXTREME  CASES 
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Case  A.  As  P -♦  «>,  the  repair  periods  have  shorter  and  shorter  durations, 
and,  in  the  limit,  repair  is  instantaneous.  Intuitively,  then,  as  P 
the  model  developed  in  Section  1 approaches  the  queue  not  subject 

to  breakdown  and  having  constant  arrival  rate  A and  service  rate  ti. 

That  this  is  in  fact  the  case  can  be  shown  from  (15),  (27),  ^nd  (51);  as 
P -^00,  G(*),  L,  and  W all  approach  the  corresponding  quantities  for 
the  H/E^l  queue. 


Case  B.  If  k = 1 and/or  m = 1,  exponential  service  and/or  repair  times 
result.  When  both  k = 1 and  m = 1,  expression  (2j)  fo  • L reduces  to 
a special  case  (4^  = O)  expression  (55)  Yechiali  and  Naor  [7,  p.  729]. 


I 


<1 


i 


Case  C.  Constant  service  and/or  repair  times  can  be  analyzed  by  letting 
k and/or  m -» <».  It  is  clear  from  (2/)  and  (5I)  that  L and  W 

are  decreasing  and  convex  functions  of  both  k and  m.  Their  limiting 
values  are  obtained  by  replacing  [(k+l)/2k]  and/or  [(m+l)/2m]  by 
1/2  in  (27)  and  (5I), 

Case  D.  If  the  Poisson  arrival  rate  is  stationary  (A^^  = A),  expressions 
(27)  and  (51)  for  L and  W reduce  to  special  cases  (Erlang  service 
and  repair)  of  relationships  {2k)  and  (26)  of  Avi-Itzhak  and  Naor 
[1,  P.  309]. 

I 
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Case  E.  In  some  practical  situations^  no  customers  enter  the  system 
when  it  is  failed^  either  by  their  own  choice  or  because  of  restrictions 
by  the  system.  During  the  repair  process,  then,  customers  neither  enter 
nor  leave  the  system.  Thus,  it  is  intuitive  that  not  only  L but  also 
the  steady-state  probabilities  of  having  n customers  in  the  system 
(n  = 0,  1,  2,  ...)  are  equal  to  those  for  the  queue  not  subject 

to  breakdown  and  having  constant  arrival  rate  A and  service  rate  4. 
That  this  is  in  fact  the  case  can  be  seen  by  setting  = 0.  Then 
r = A/|j  and,  from  ( I5)  and  (27),  both  G{z)  and  L simplify  to  the 
corresponding  quantities  for  the  M/E^/l  queue.  Of  course,  as  (Jl)  with 
Aj^  = 0 indicates,  W is  greater  than  in  the  case. 


, 3.  BEHAVIOR  OF  L AND  W AS  A FUNCTION  OF  c 

! In  order  to  investigate  how  the  system  behaves  as  a function  of  c, 

I the  assumption  of  Sections  1,  2,  and  1+  that  c = 1 is  now  dropped.  Recall 

that  c controls  how  rapidly  the  system  occilates  between  the  operational 
and  failed  states.  Varying  c while  holding  0!  and  p constant  does 

j not  change  p and  p the  steady-state  probabilities  of  the  system 

I 

I being  operational  and  failed,  respectively.  However,  as  c increases, 

the  system  fluctuates  more  rapidly  between  the  operational  and  failed 
states,  or,  equivalently,  the  mean  time  the  system  stays  in  each  state 
I approaches  0. 

4 

I 

( 

I 

I 
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Expressions  for  L and  tf  as  a function  of  c can  be  obtained 
by  replacing  cc  and  p by  cCC  and  cp  everywhere  in  (2j)  and  (51) , ^ 

From  the  expressions  that  result,  it  is  easily  shown  that 

L'(c)  = -ac 

-2  ' 

W'(c)  = -be 

I 

where  a and  b are  positive  constants  involving  A,  Tv,,  ^ a and  p, 

I 

Hence,  both  L and  W are  decreasing  and  convex  functions  of  c.  ' 

To  interpret  this  result  qualitatively,  consider  two  equally 
reliable  systems  (i.e,,  identical  p 's)  also  having  identical  k,  m.  A,  I 

i 

and  however,  suppose  one  system  has  infrequent  failures  but  long 
repair  times  (a  low  c)  and  the  other  undergoes  frequent  but  quickly 
repaired  failues  (a  high  c) . If  the  objective  is  to  minimize  L and  W, 
then  the  latter  system  should  be  chosen. 

On  a more  quantitative  level,  the  result  supports  a general  con- 
jecture of  Ross  [6]  that,  in  a single  server  infinite  capacity  queueing 
model,  the  "more  stationary"  the  Poisson  arrival  process  is,  then  the 
smaller  the  average  customer  delay.  As  in  [2]  and  [6],  this  conjecture 
has  been  verified  in  a special  case.  To  see  this,  the  behavior  of  the 
system  as  c ->  « will  be  investigated.  It  is  easy  to  show  from 
(27)  and  (31)  that  as  c -» »,  both  L and  W approach  the  corresponding 
quantities  for  an  queue  not  subject  to  breakdown  and  having  con- 

stant  arrival  rate  A and  service  rate  |i.  Thus  as  c -»oo,  the  system 


becomes  more  stationary  in  the  sense  that  it  behaves  more  and  more  like 


the  m/E  /l  queue  with  parameters  A and  Also  note  that  since  L 

k 


and  W are  decreasing  in  c,  the  smallest  values  they  can  ever  achieve 
are  the  corresponding  values  for  the  m/Ej^^I  queue  with  parameters  A 

/V 

and 


APPLICATION  TO  A PRODUCTION- STORAGE  SYSTEM 


By  regarding  the  server  as  a production  process  turning  out  items 
one  at  a time  and  each  customer  as  a unit  demand  for  the  product^  the 
results  of  Section  1 can  be  used  to  analyze  a production- storage  process 
subject  to  breakdown  and  having  the  following  additional  characteristics: 

(a)  Unsatisfied  demand  is  always  backlogged. 

(b)  Items  not  needed  immediately  to  satisfy  backlogged  demand 
are  stored  for  future  use  up  to  a level  of  S^  the  finite  capacity  of 
the  storage  facility. 

(c)  When  the  storage  facility  is  filled,  no  production  takes  place. 


1 


li 


I ’ 


In  such  a production- storage  model,  two  quantities  of  interest  are  I, 
the  average  inventory  level  where  a negative  value  indicates  that  the 
items  are  backlogged  rather  than  in  storage,  and  R,  the  fraction  of  time 
that  demand  can  be  met  without  backlogging.  For  example,  it  may  be  desired 
to  choose  S so  that  I and/or  R are  greater  than  some  specified 
design  parameters.  Meyer  Rothkopf  and  Eldrcdge  [5]  consider  four 


production- storage  models  differing  from  the  one  just  described  in  that 
production  and  demand  both  occur  at  constant  rates,  unsatisfied  demand 
is  lost,  and  the  durations  of  operating  and  failed  periods  have  distri- 
butions corresponding  to  the  four  possible  combinations  of  constant  and 
exponential. 

A state  (i,j)  of  the  queueing  system  translates  into  a state 
of  the  production- storage  system  as  follows:  0 < j < kS  is  equivalent 
to  kS-j  phases  in  storage  and  j > kS  is  equivalent  to  j-kS  phases 
backlogged.  Of  course,  i has  the  same  interpretation  in  both  models. 
Given  this  one-to-one  correspondence  between  states  in  the  two  models, 
expression  {2rf)  and  recursions  ( l6)  can  be  used  to  compute 


I = S-L 


R = V 


.m 


vk(s-i) 


^i=0  ^j=0 


Pij 


(S  > 0) 


where  a negative  value  for  I indicates  a backlog.  As  a numerical 
example,  consider  a system  with  a = 1,  P = 5,  m = 5,  X = 8, 

= ^,  4 = 16,  , k = and  suppose  it  is  desired  to  choose  S so  that 

both  I >5  and  R >0.95«  Computed  with  the  aid  of  (?7)  and  (16),  the 
following  table  shows  I and  R as  a function  of  S: 


i-i 


I 


I 

I 


i 


I 
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APPENDIX: 


DERIVATION  OF 


H(y) 


To  obtain  H(y)  multiply  each  equation  of  the  sets  {(la)^  (Id)j^ 
{(lb),  (le)),  and  {(Ic),  (If))  by  y"  where 


n = 


0 

k 

j modulo  k 


if  j = 0 

if  j k,  2k,  3k,  . . . 
otherwise  . 


Summing  over  all  j then  yields,  respectively. 


Hjy)  = (mp)'^  [aHo(y)  - ^ (52) 

\ PiO^^-y*'^  (1  < i < ">-!)  (55) 

HQ(y)  = [ay  - kn(  1-y)  ]‘^  {m^y  Hj^(y)  - knPpQ(  1-y)  - y(  l-y'^) 

' ^^Poo  ^ - Pol^^^  • (5^) 


Equations  (32)  and  (33)  c^n  be  used  recursively  to  express  H^(y)  in 
terms  of  ^^(y)  as 


Hj(y)  = (inP)’^  [QHo(y)  - Ai(i-y'^)  7™^^ 


(35) 


K 


■ f 

and  (^h)  can  be  rearranged  as 

Hj^(y)  = (mPy)'^  I [ay  - k4(  l-y)  ] H^(y')  + k^i  Poq(  1-y) 

+ y(i-y*')  [?^Pqq  + 

Equating  the  expressions  for  Hj^(y)  in  (35)  ^nd  (36)  solving  for  ^^(y) 
results  in 

^o^y)  = pQo  [[y(i-y'^)  V[Mi-y)]) 

{Ai  Pj^Q  + APqq  + k4[HQ(0)  - Pqj^])  . (37) 

Using  (32)  and  (33)  recursively  yields 

H^{y)  = (^)'^  CQ«o(y)  - A^(l-y‘^)  Pno^  . (1  < i <">)  (58) 

H(y)  can  now  be  obtained  from  (24),  (37),  (38)  as 

H(y)  . Kc^mi  { Poo  - life}  '^"oo  ♦ \ J'°=1  "lo  • - "oi)  >} 

- (A^/ii<i)  (l-y*")  E"^^  ip.Q  . (39) 

Using  H(  1)  = 1 and  computing  lim^  H(y)  by  applying  L'Rdpital's 
rule  to  (39)  results,  after  algebraic  simplification^  in 
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kn[Ho(0)  - P.Q  - A Pqq  . (W) 


Substituting  (UO)  into  (39)  gives 


H(y)  = [(CtfP)/p]pQQ  + r{[y(l-y  )]/[k(l-y)]) 


- (X^/r<^)  (l-y^)  ip.Q 


(Al) 


Equations  ( la)  and  ( lb)  can  be  rewritten  as 


PnO  = ""Poo  - ^1  PtrO 


^ PiO  = Pi+1  0 ■ ^1  PiO  ( 1 < i < m-l) 


which  is  equivalent  (by  induction)  to 


PiO  = (%0  • h Ci  PnO^  ">-1,  1).  (^2) 


Using  (42)  to  evaluate  Z™_q  P^q  yields  the  relationship 


1-p  = [(Q^t3)/p]pQQ  - (Aj^/mP)  F”^]^  iPj^Q  . 


(i^3) 


Because  of  (43)  and  (14),  (41)  simplifies  to  (25)^  the  expression  for 
H(y)  given  in  Section  1. 
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